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摘 要 : 针对 生态 环境 脆弱 的 寒 旱 区 


城乡 建设 .生态 环境 监测 与 保护 等 具有 重要 意义 


展 地 物 要 素 提 取 以 及 土地 覆盖 变化 监测 研究 ,对 农业 规划 、 
o TEE} 2015—2019 4 3r iE 4 2 S E X Landsat-8 


影像 构建 数据 集 , 对 比 3 种 传统 方法 :最 大 似 然 分 类 (Maximum likelihood classification, MLC) 支持 
向 量 机 (Support vector machine,SVM) 和 随机 森林 (Random forest,RF) 及 5 种 语义 分 割 模型 Deep- 
Labv3 + (Xception ) , DeepLabv3 + (MobileNet ) .SegNet (ResNet50) , U- Net (MobileNet ) 和 PSPNet (Mo- 
bileNet) ,选取 最 优 自动 化 地 物 提取 模型 对 研究 区 1998 一 2020 年 农用 地 、 建 筑 用 地 、 水 体 和 荒漠 4 种 


地 物 要 素 进行 


分 类 ,并 运用 土地 利用 转移 矩阵 和 动态 度 进行 定量 动态 变化 分 析 。 结 果 表 明 :Deep- 


Labv3+(Xception ) 模 型 可 以 实现 更 准确 .更 高 效 的 地 物 提取 ,总 体 精 确 度 (OA) Kappa 系数 和 Fl 值 
分 别 为 96.06% 、0.96 和 0.86, 其 中 所 选 模型 的 平均 交 并 比 (MIoU ) 较 其 他 模型 提升 0.03~0.39。 近 23 


a 莫 索 湾 垦区 的 荒漠 、 农 用 地 和 建筑 用 地 三 者 的 土地 结构 转化 较为 明显 ,荒漠 总 
积 增加 2.33% ,水 体 


农用 地 总 面积 增加 12.68% ,建筑 用 地 总 


Hye > 15.00% ， 
积 变 化 较为 平稳 。 地 物 类 型 总 体 


转变 方向 为 荒漠 向 农用 地 转化 、 农 用 地 向 建筑 用 地 转化 。 该 研究 可 为 深度 学 习 技术 应 用 于 中 分 辩 
率 遥 感 卫星 影像 领域 中 实现 土地 利用 及 变化 动态 监测 提供 参考 。 
关 键 词 : 遥感 影像 ;长 时 间 序 列 ; READ; DeepLabv3+; 地 物 分 类 ; 动态 监测 


文章 编号 : 


遥感 技术 的 迅速 发 展 为 土地 利用 及 履 盖 变化 
(Land use and cover change, LUCC) #é( T ria", 
对 生态 环境 脆弱 的 干旱 区 土地 利用 及 其 变化 情况 
准确 及 时 监测 提供 重要 支持 ”。 新 疆 为 典型 干旱 
区 , 莫 索 湾 暴 区 地 处 新 疆 天 山北 坡 经 济 带 的 重要 发 
展区 ,为 荒漠 区 绿洲 型 农作物 区 , 近 些 年 随 着 人 类 
社会 活动 加 出 ,绿洲 内 外 面临 还 漠 化 等 生态 问题 。 
因此 ,快速 准确 实现 该 地 区 土地 覆盖 及 变化 分 析 进 
而 获取 空间 分 布 数据 ,为 环境 脆弱 区 域 的 生态 环境 
建设 5 城乡 规划 和 农业 管理 等 提供 有 效 技术 
手段 。 

遥感 影像 分 类 方法 中 ,基于 像 元 分 类 的 传统 机 
器 学 习 算 法 ,如 最 大 似 然 分 类 (Maximum likelihood 
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classification , MLC ) 、 文 持 向 量 机 (Support vector ma- 
chine, SVM) 和 随机 森林 (Random forest, RF) 等 , 借 
助 影像 纹理 .光谱 等 信息 进行 分 类 , 易 产 生 “ 同 物 异 
谱 "”“ 同 谱 异 物 " 以 及 椒盐 噪声 现象 "”。 面 向 对 象 
影像 分 析 方 法 (Object-based image analysis, OBIA ) 
面临 人 为 参与 分 割 参数 调整 ,缺乏 上 下 文 指导 等 难 
题 *“"。 近 年 来 ,深度 学 习 技 术 在 挖掘 海量 遥感 数 
据 方面 的 潜力 逐渐 凸显 "中 。2015 年 ,Long 等 中 提 
出 全 卷 积 神经 网 络 (Fully convolutional networks, 
FCN) ,将 全 连接 层 蔡 换 为 卷 积 层 实现 更 优 特征 提取 ， 
后 续 出 现 诸多 基于 FCN 改进 的 语义 分 割 模型 "1 分 
制 和 识别 精度 不 断 提高 。 

DeepLab 系列 ”模型 在 不 断 更 新 ,最 新 的 
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DeepLabv3+ 模 型 采用 编 解码 架构 和 空洞 空间 金字 
塔 池 化 模块 (Atrous spatial pyramid pooling, ASPP) , 
能 够 结合 上 下 文 信息 提取 多 尺度 特征 图 ,不 断 恢 复 
空间 信息 ,克服 了 传统 方法 特征 表达 不 全 缺乏 上 
下 文 指导 等 不 足 。 学 者 们 借助 语义 分 割 网 络 模型 ， 
完成 了 单一 或 多 地 物 要 素 分 类 的 研究 ””。 许 慧 敏 
等 中 借助 U-Net 模 型 的 压缩 通道 和 扩展 通道 达到 较 
高 定位 精度 和 识别 精度 ,优化 了 分 类 效果 。 杨 建 宇 
等 2 通过 SegNet 融 合 编 解码 (Encoder-decoder) 架 构 
及 跳跃 网 络 ,在 提取 农村 建筑 物 用 地 中 Fl 值 高 达 
0.91。 刘 文雅 等 ”实验 表明 DeepLabv3+ 长 于 提取 高 
层 特 征 , 获 得 总 精度 91.02% 最 优 城市 绿地 要 素 提取 
效果 。 由 此 可 见 深度 学 习 能 够 充分 挖掘 遥感 影像 
中 的 光谱 及 空间 信息 ,在 实现 影像 自动 化 特征 提 
取 , 快 速 准确 获取 地 物 信息 方面 有 较 大 潜力 。 
在 土地 利用 变化 监测 方面 ,学 者 运用 传统 机 带 
学 习 算 法 的 分 类 结果 实现 长 时 间 序 列 的 变化 监测 ， 
由 于 传统 机 顺 学 习 算 法 自动 化 程度 低 ,在 长 时 间 序 
列 影 像 分 类 任务 中 人 机 交互 现象 频繁 ,分 类 精度 和 
言 息 提 取 效 率 都 有 进一步 提升 空间 ””。 另 有 学 者 
不 断 尝 试 将 深度 学 习 技 术 应 用 于 遥感 影像 监测 土 
地 利用 变化 方面 ,多 数 是 借助 卷 积 神 经 网 络 实现 双 
时 相 变 化 检测 3, 但 利用 深度 学 习 实 现 多 时 相 匠 
感 影 像 地 物 分 类 后 变化 监测 的 研究 仍 未 深入 开 
展 。Liu 等 ”借助 FCNs 获取 印度 班加罗尔 2012 一 
2016 年 贫民 窒 土 地 稚 盖 信息 后 分 析 动 态 监 测 平民 
富 变 化 状况 ,但 该 变化 分 析 研 究 的 影像 数据 周期 
综 上 ,目前 针对 莫 索 湾 垦 区 土地 利用 结构 及 其 
动态 变化 检测 还 有 待 进一步 利用 新 工具 和 方法 深 
入 开展 研究 。 人 工 进 行 实地 地 物 监测 成 本 高 效率 
低 , 并 存在 信息 更 新 渍 后 性 ,运用 传统 算法 自动 化 
程度 不 高 ,影响 地 物 信 息 获 取 的 效率 与 准确 度 ,而 
将 遥感 影像 数据 与 深度 学 习 算 法 结合 实现 自动 化 
特征 提取 更 为 有 效 可 行 ,并 且 跨 时 间 长 度 超过 20 a 
的 相关 研究 偏 少 ”。 因 此 ,本 文采 用 Landsat 系列 
(CTM,ETM+,OLIU 卫 星 影像 作为 源 数据 ,制作 地 物 分 
类 数据 集 ,对比 训练 3 种 传统 方法 (MLC .SVM 和 
RF) 和 5 种 语义 分 割 模型 [DeepLabv3+(Xception)、 
DeepLabv3 + (MobileNet) , SegNet (ResNet50) , U-Net 
(MobileNet) 和 PSPNet(MobileNet) ] ,选择 表现 最 佳 
地 物 类 型 提取 网 络 实现 1998 一 2020 年 多 时 相 地 物 


分 类 ,并 进一步 分 析 地 物 结构 变化 及 其 驱动 力 
素 ,为 莫 索 湾 明 区 为 例 的 寒 旱 区 自然 资源 时 空 变化 
分 析 和 政府 宏观 决策 提供 参考 。 


1 研究 区 与 数据 来 源 


1.1 研究 区 概况 

Wt FE RK BER YS BK (44° 23'~45° 12’N, 85°52'~ 
86° 19'E) Hb ch ir He Ae LCRA PE , TA AS DE TA BE 
沙漠 南 缘 , 沿 主干 道 石 莫 公 路 分 布 4 个 团 场 :147 
十 户 滩 镇 .148 团 西 营 镇 149 团 东 旱 城 镇 和 150 
西 古城 镇 及 其 周边 附属 (图 1) ,实地 覆盖 面积 约 
14500 km?。 其 中 147 团 主要 分 布 有 优先 发 展 的 化 
工 新 材料 工业 区 ,148 团 .149 团 和 150 团 农业 发 展 
结构 以 棉花 种 植 为 主 ,林果业 和 设施 农业 为 辅 。 研 
究 区 属于 典型 温带 大 陆 性 干旱 气候 ,地 势 平 坦 , 作 
为 灌溉 农业 区 ,其 耕地 资源 丰富 ,是 新 疆 棉 花 的 主 
产 区 之 一 ,城镇 分 布 较为 密集 ,产业 结构 合理 ,为 重 
BAR RKO, 


图 1 rs TLR BX 


Fig. 1 Map of Mosuowan reclamation area in Xinjiang 


1.2 实验 数据 
实验 数据 包括 Landsat 系列 卫星 影像 数据 ( 表 
1) 及 莫 索 湾 明 区 的 土地 利用 矢量 数据 ,其 中 卫星 影 
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不 盼 丽 等 ,基于 深度 学 习 的 寒 旱 区 多 时 序 影像 土地 利用 及 变化 监测 一 一 以 新 疆 莫 索 湾 晨 区 为 例 


表 1 1998—2020 Ẹ Landsat 系列 卫星 影像 数据 详情 
Tab.1 Details of Landsat satellite image data from 1998 to 2020 


传感器 类 型 发 射 时 间 (年 -月 -日 ) 空间 分 辨 率 /m 时 间 分 辨 率 /4 pane 云 量 范围 /% 
Landsat-5/TM 1984-03-01 30 16 1998 ,2006—2011 0.00~3.00 
Landsat-—7/ETM+ 1999-04-15 30 16 1999—2005 0.00~6.45 
Landsat-8/OLI 2013-02-11 30 16 2013—2020 0.00~2.39 
像 数 据 取 自 地 理 空 间 数 据 云 官网 (https://www.gs- 开始 
cloud.cn) 及 美国 地 质 勘 测 局 (https://earthexploere. = ,-------------------Y---------------- 
usgs.gov) 官 网 。 卫 星 重 访 周期 短 (16 d) 增 大 了 实验 Landsat-5 || Landsat~7 || Landsat-s 
数据 的 可 选择 性 。 研 究 区 作物 生长 季 (7 一 10 月 ) 影 数 
被 光谱 和 纹理 特征 易 与 其 余 ATE KAY 据 
像 植被 光谱 和 纹理 特征 易 与 其 余地 物 特征 区 分 i mN 
选择 云层 覆盖 量 低 于 3%(2010 年 一 景 数据 云 量 达 z 等 前 期 处 理 
6.45% ) 的 生长 季 影 像 数据 ,选取 的 多 时 相 过 感 影像 M [a ER 
质量 较 高 。 整 幅 影 像 大 小 为 1133x2874 像素 ,数据 
paa a 
时 间 跨 度 23 a, 共 22 景 。 
a 训练 数据 测试 数据 
"E o ee 
nooo (expla 
fir IEH 示 。 ,对 研究 区 Landsat SegNet、U-Net 等 ) 
技术 流程 如 图 2 所 首先 ,对 研究 区 Land 
系列 卫星 影像 进行 预 处 理 并 制作 数据 集 ,开展 传统 a 
方法 和 语义 分 割 模型 对 比 实验 后 ,选取 分 类 效果 最 N 传统 方法 
好 的 模型 ,以 保证 土地 分 类 提取 结果 的 可 靠 性 ;第 Bt ‘sci eit 
二 步 , 借 助 地 物 类 型 提取 模型 完成 1998 一 2020 年 莫 
步 借助 也 物 类 提取 模 完成 a O] 待 分 类 数据 
索 湾 星 区 4 种 典型 地 物 要 素 的 分 类 ,进一步 验证 本 7 7 
文 方法 的 有 效 性 及 适用 性 ,并 绘制 地 物 分 类 专题 图 地 物 类 型 提取 模型 
可 视 化 分 类 效果 ;最 后 ,根据 分 类 结果 计算 土地 利 Fr---------- 二 -一 -一 一- 一- 一- 一- 二 -一 - 
用 转移 矩阵 和 动态 度 ,分析 以 莫 索 湾 明 区 为 例 的 寒 [地 物 类 型 提取 结果 待 分 类 数据 集 | 
= 23 Tal Ae a Of: HE Fy ZB e We al 类 y 
旱 区 各 地 物 要 素 空 间 分 布 特点 并 进行 变化 监测 2 
分 析 。 T CE y y 
22 数据 集 创 建 及 样本 组 织 方式 上 [oalfprecision}recatl} [F1] | MioU] | Kappa 
| 
2.2.1 预 处 理 ” 数 据 预 处 理 步 又 分 为 :辐射 定 
。 数据 页 处 数据 预 处 akin 辆 射 定 A | 
标 、 大 气 校 正 .研究 区 裁剪 .数据 格式 转换 等 。 其 中 
2003 一 2005 年 3 景 Landsat-7 ETM+ 影 像 数据 出 现 条 o| 土地 转移 矩阵 、 动 态度 “|>| 变化 监测 分 析 


带 数据 丢失 ,借助 ENVI 5.5.2 软 件 扩展 工具 land- 
sat_gapfill 完成 像素 空间 插值 修补 ,条 带 修复 后 的 数 
据 几乎 不 影响 分 类 效果 。 

用 于 模型 训练 和 测试 的 数据 为 2015 一 2019 年 
间 5 景 影像 数据 ,由 于 整 幅 影像 尺寸 过 大 , 受 实验 软 
硬件 性 能 限制 ,输入 语义 分 割 模 型 的 数据 需 先 进行 


注 :OA 为 总 体 精确 度 ;precision 为 查 准 率 ;recall 为 查 全 率 ;F1 为 F1 
值 ;MIoU 为 平均 交 并 比 ;Kappa 为 Kappa 系 数 ;SVM 为 支持 向 量 
机 ;RF 为 随机 森林 ;MLC 为 最 大 似 然 方 法 。 

图 2 地 物 分 类 及 覆盖 变化 动态 监测 技术 流程 图 


Fig. 2 Flow chart of dynamic monitoring technology for 


feature classification and coverage change 
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PEAR) BAD EE ,综合 考虑 切割 大 小 定 为 416x416 像 
素 , 滤 除 空白 影像 块 ,一 张 影像 切割 成 20 张 相互 独 
立 影像 块 ,切割 完成 后 共 100 个 影像 块 。 

2.2.2 数据 集 制 作 鉴于 Landsat 系 列 影 像 数据 30 
m 空 间 分 辨 率 以 及 研究 区 地 物 分 布 情况 ,参考 《土地 
利用 现状 分 类 》(CB/T21010-2017) 标 准 , 将 研究 区 
的 地 物 类 型 分 为 5 类 ,分 别 是 建筑 用 地 、 农 用 地 、 水 
体 .荒漠 以 及 背景 ( 表 2)。 综 合 分 析 影 像 的 光谱 信 
息 和 纹理 特征 总 结 了 以 下 解 译 标 志 : 城 乡 建筑 群 呈 
现 紧 凑 块 状 , 影 像 大 多 呈 和 白色 ,部 分 影像 有 红色 E 
色 夹 杂 , 呈 点 状 散布 ;荒漠 的 影像 呈 褐 色 .黄色 , 呈 
不 规则 片 状 ;农用 地 影像 为 墨绿 色 , 且 多 呈 大 片 连 
续 不 规则 状 ; 水 体 颜色 呈 深 浅 不 一 的 绿色 ,分 布 区 
域 较为 固定 ,为 连续 片 状 , 边 界 较为 明显 。 结 合 高 


表 2 研究 区 地 物 类 型 信息 
Tab.2 Land use type information in the study area 


颜色 编码 


地 物 类 型 包含 类 型 (Red , Green, Blue) 
建筑 用 地 ”城镇 ,村 庄 .工业 建筑 用 地 等 (1,1,1) 
农用 地 耕地 (旱地 水 田 、 光 果园 等 )、 林 (2,2,2) 
地 (灌木 林地 、 有 林地 等 ) .草地 
(天 然 牧草 地 .其 他 草地 ) 等 
水 体 滩涂 . 坑 塘 .水库 ,湿地 等 (3,3,3) 
荒漠 沙漠 未 利用 土地 等 (4,4,4) 
背景 研究 区 范围 以 外 的 区 域 (0,0,0) 
编码 部 分 [了 
1x1 卷 积 DA 
深度 卷 积 网 络 || ZSEB |> 


模型 输入 


| 
| 
| 
| 
H 
| 
| 
| 
| 
| 


低层 特征 


| 

| 4 

| p -9- 份 3x3 卷 积 >| 四 倍 上 采样 
| 

| 


4) HEE Google Earth 影像 进行 目 视 解 译 , 在 LabelMe 
深度 学 习 标 注 工 具 中 完成 2.2.1 节 中 影像 块 的 数据 
标注 工作 ,进而 转化 为 颜色 编码 (Red, Green, Blue) 
的 数据 集 。 

2.2.3 数据 样本 组 织 方式 ”为 解决 样本 不 足 带 来 过 
拟 合 问题 ,利用 数据 增强 方式 将 2.2.2 节 中 人 工 标注 
的 数据 集 进 行 扩充 ,以 提高 样本 多 样 性 。 第 一 阶段 
数据 增强 在 模型 训练 前 采用 图 像 30° 60° 90° .180” 
旋转 、 高 斯 噪声 扰动 以 及 小 样本 复制 方式 ,将 100 张 
数据 扩充 到 630 张 ;第 二 阶段 数据 增强 是 模型 训练 
过 程 中 ,在 数据 生成 器 里 通过 色彩 抖动 .比例 缩放 、 
左右 翻转 方式 随机 增强 数据 ,进一步 增强 了 模型 鲁 
棒 性 。 最 终 用 于 模型 训练 的 数据 共 306 张 ,用 于 模 
型 验证 的 数据 共 204 张 ,用 于 模型 测试 的 数据 共 
120 张 。 

2.3 DeepLabv3+t =! 

本 文 使 用 的 DeepLabv3+ 架 构 如 图 3 所 示 ,该 模 
型 是 由 Chen 等 2 提出 的 编码 -解码 架构 。Deep- 
Labv3+ 模 型 作为 编码 器 ,采用 先 串 行 空 洞 卷 积 后 并 
行 空洞 卷 积 架 构 进行 特征 提取 ,不同 扩 张 率 并 行 空 
洞 卷 积 可 实现 不 同 尺 度 特征 重 采 样 进 而 编码 全 局 
上 下 文 特征 ,空洞 卷 积 通过 扩张 率 增 大 感受 野 (Re- 
ceptive field ) 的 同时 不 增加 参数 量 ,实现 精度 与 计算 
速度 的 平衡 。 新 增 解 码 器 部 分 ,编码 器 输出 的 特征 


一 :一 一 :一 :一 :一 :一 :一 一: 一: 


| 
| 
| 
| 
| 
æ 
| 
| 
| 
| 


一 一 
| 
| 


图 3 DeepLabv3+ 模 型 架构 图 
Fig. 3 DeepLabv3+ model architecture 
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图 作为 该 部 分 的 输入 ,采用 4 倍 双 线 性 插值 上 采样 ， 
得 到 的 特征 信息 与 串 行 空 洞 卷 积 的 输出 进行 向 量 
拼接 ,实现 高 层 与 低层 特征 融合 ,再 经 过 常规 3x3 卷 
积 和 4 倍 双 线性 插值 上 采样 后 得 到 分 类 结果 。Dee- 
pLabv3+ 模 型 将 空洞 卷 积 与 深度 分 离 卷 积 结合 成 空 
洞 可 分 离 卷 积 ,应 用 于 编码 器 的 ASPP 模 块 和 解码 
结构 中 ,二 者 结合 保证 分 类 精度 的 同时 减少 运 
算 量 。 
2.4 实验 环境 及 模型 参数 配置 

本 文 涉及 的 5 个 网 络 模 型 均 在 Keras (后 端 为 
Tensorflow ) 深 度 学 习 平 台 下 搭建 ,计算 机 硬件 平台 
为 16 核 Inter(R)Xeon(R )Sliver 4110 CPU@2.10GHz, 
GPU 为 NVIDIA Quadro P4000,32G 内 存 。 语 义 分 割 
网 络 模型 主要 参数 设置 :学 习 率 (Learning rate) 控 制 
模型 学 习 速率 ,为 0.0001; 学 习 率 变化 指数 (Camma ) 
控制 学 习 率 变化 速率 ,为 0.5; 人 迭代 次 数 (Epoch ) 为 
500; 24.52 E (Patience ) 为 8, 当 学 习 率 有 8 次 迭代 不 
变 时 学 习 率 更 新 ; 批 处 理 尺寸 0Batchsize ) 为 3; 优 化 
fire HY Adam. 
2.5 模型 精度 评价 指标 

选取 基于 分 类 混淆 和 矩阵 (Confusion matrix ) 运 算 
的 6 种 模型 精度 评价 指标 ,分 别 为 总 体 精 确 度 
(Overall accuracy, OA) 、 查 准 率 (precision,%)、 查 全 
率 (recall,%)、Kappa 系数、 平均 交 并 比 (Mean Inter- 
section-over Union, MIoU) fil F1 值 (Fl-score) ,其 中 
Kappa 系 数 用 于 检验 模型 预测 结果 与 真实 值 一致 性 
的 指标 , MIoU 是 基于 全 局 评价 的 ,Fl 值 为 查 全 率 和 
查 准 率 的 调和 平均 数 。 上 具体 计算 公式 如 下 : 


Tp TN +P FEN a 

precision = aF (2) 

recall = PTN (3) 

poe aE (4) 
mxm 

Kappa= 让 二 E (5) 

F1 - score = 2 (6) 


1 1 


ean 
precision recall 


_ TP 
MloU = spi (7) 


式 中 :TP(True positive ) X FLIE Pi]; TN (True negative ) 
HAH ill ; FP (False positive) 为 假 正 例 ;FN (False 


negative) 为 假 负 例 ; P, 为 每 一 类 正确 分 类 的 样本 数 
量 之 和 除 以 总 样本 数 , 即 总 体 分 类 精度 ; a, 为 第 7 
类 的 真实 样本 个 数 ; b, 为 预测 出 来 的 第 nn 类 样本 个 
Bl; n 为 共有 n 类 地 物 类 型 ; m 为 所 有 地 物 分 类 中 
总 共有 m 个 像 元 数 。 

2.6 土地 利用 及 覆盖 变化 分 析 

2.6.1 Let] ABBR 土地 利用 转移 算 阵 反应 
不 同时 段 内 地 物 类 型 相互 转换 关系 ,揭示 不 同 地 物 
类 型 间 的 转移 速率 ,是 马尔 可 夫 模 型 在 土地 利用 及 
变化 监测 方面 的 应 用 。 借 助 影像 分 类 结果 运用 EN- 
VI 软件 中 的 土地 利用 类 型 转移 矩阵 ,对 1998 年 和 
2020 年 2 期 分 类 影像 完成 土地 利用 转移 统计 ,土地 
利用 转移 矩阵 公式 如 下 : 


CI Ay a, 
fala “Se ay (8) 
Q mi a m2 S Q mn 


式 中 : 4 为 转移 和 矩阵; a,, 为 7 时 期 的 m 种 地 物 类 
型 转变 为 T+1 时 期 n 种 地 物 类 型 的 面积 (km ); m 、 
n 分 别 为 初始 研究 期 和 结束 研究 期 的 地 物 类 型 ,其 
中 为 2 个 时 期 地 物 类 型 没有 发 生 转 变 。 

2.6.2 动态 度 ”借助 土地 利用 动态 度 定 量 展示 人 研究 
区 的 地 物 类 型 动态 演变 过 程 .变化 速度 以 及 变化 程 
度 汪 。 单 一 土地 动态 度 公式 如 下 : 


OO 
~ U 


K x A x 100% (9) 


式 中 : K 为 特定 时 间 范 围 内 某 一 地 物 的 动态 度 , 即 
年 变化 率 (km…a'); U,、U 分别 为 研究 初 斯 和 研究 
末期 某 一 地 物 总 面积 (km); T 为 变化 时 间 段 (a)。 


3 结果 与 分 析 


3.1 地 物 要 素 提 取 分 类 实验 

3.1.1 模型 分 类 结果 对 比 ” 为 证 明 所 选 模型 的 有 效 
性 ,在 相同 测试 集 上 对 3 种 传统 算法 和 5 个 语义 分 
割 模型 进行 分 类 结果 对 比 ( 表 3)。 由 于 传统 机 带 学 
习 方 法 忽略 影像 深层 语义 信息 , 且 依 赖 于 人 机 交 
互 , 因 此 总 体 上 语义 分 割 模型 在 分 类 精度 上 优 于 传 
统 机 器 学 习 算 法 。 而 深度 学 习 技 术 通 过 输入 数据 
集 对 特定 分 类 模型 进行 训练 ,借助 训练 好 的 模型 实 
现 大 批量 遥感 影像 分 类 工作 ,其 自动 化 程度 高 于 传 
统 方法 ,并 且 能 更 好 地 挖掘 高 层 语义 信息 ,分 类 效 
果 更 优 。DeepLabv3+(Xception) 模 型 对 地 物 分 类 有 
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3 不 同 模型 各 项 评价 指标 


Tab. 3 Evaluation indexes of different models 


ae 评价 指标 

OA/% Kappa precision/% recall/% F1 值 MIoU 
MLC 62.73 0.42 52.78 71.17 0.51 0.38 
SVM 68.32 0.48 52.12 73.18 0.54 0.41 
RF 76.35 0.59 75.59 85.69 0.75 0.64 
DeepLabv3+ (Xception) 96.06 0.96 87.69 83.78 0.86 0.77 
DeepLabv3+ (MobileNet) 94.77 0.94 74.59 75.82 0.75 0.62 
SegNet (ResNet50) 95.47 0.95 81.94 81.04 0.81 0.71 
U-Net (MobileNet) 95.83 0.95 85.23 83.13 0.83 0.74 
PSPNet (MobileNet) 95.52 0.90 85.34 81.08 0.83 0.74 


RF 为 随机 森林 。 下 同 。 


明显 的 精度 优势 ,分 类 后 的 总 体 精 确 度 Kappa 系 
数 、 查 准 率 、 查 全 率 、Fl 值 .MIoU 分 别 为 96.06% 、 
0.96 ,87.69% 、83.78% 、0.86、0.77, MIoU 比 其 他 模型 
提升 了 0.03~0.39 ,表明 该 网 络 在 中 分 辩 率 卫星 遥感 
影像 中 对 典型 地 物 要 素 分 类 提取 能 力 优 于 其 他 模 
型 ,适用 性 较 好 。 

表 4 展 示 了 各 个 网 络 上 4 种 地 物 要 素 提 取 对 比 
结果 。 总 体 上 ,本 文 方法 对 中 分 辩 率 遥感 影像 地 物 
要 素 提取 表现 最 优 ,尤其 是 在 荒漠 地 物 形 状 不 定 ， 
建筑 用 地 呈现 出 小 面积 多 分 散 形态 下 ,各 个 模型 上 
荒漠 和 建筑 用 地 的 查 准 率 ` 查 全 率 和 Fl 值 均 低 于 
DeepLabv3+(Xception ) 模 型 ,DeepLabv3+(Xception ) 
模型 识别 建筑 用 地 的 查 准 率 、 查 全 率 、F1 值 分 别 为 
81.12% 、70.36% .0.74, 凸 显 其 在 不 同 尺 度 地 物 分 割 


注 :OA 为 总 体 精 确 度 ;Kappa 为 Kappa 系数 ;precision 为 查 准 率 ;recall 为 查 全 率 ;MIoU 为 平均 交 并 比 ;MLC 为 最 大 似 然 法 ;SVM 为 支持 向 量 机 ; 


方面 的 优势 。 

为 验证 DeepLabv3 + (Xception ) 的 适用 性 ,选取 
非 数据 集 时 期 的 影像 , 即 1998 年 和 2020 年 影像 进 
行 目 视 解 译 数据 标注 ,将 人 工 标注 的 数据 与 模型 预 
测 结果 计算 混淆 矩阵 完成 模型 精度 验证 ( 表 5)。 验 
证 结果 表明 ,1998 年 和 2020 年 的 验证 测试 数据 接近 
模型 实验 数据 ,证 明 本 文选 取 的 模型 具有 一 定 的 适 
用 性 。 
3.1.2 帝 感 影像 上 分 类 效果 对 比 5 种 模型 在 相同 
数据 集 上 的 分 类 效果 如 图 4 所 示 , 展 示 了 研究 区 局 
部 原始 影像 .标签 数据 以 及 5 种 网 络 模型 的 分 类 结 
果 。 根 据 结 果 图 对 比 ,由 于 DeepLabv3+(Xception) 
网 络 采 用 多 尺度 特征 融合 策略 ,增强 不 同 尺度 地 物 
要 素 的 提取 效果 ,对 边界 特征 的 划分 清晰 完整 , 善 


表 4 不 同 分 类 模型 对 各 要 素 分 类 结果 对 比 


Tab. 4 Comparison of classification results of different classification models for each element 


e 语义 分 制 模型 
地 物 类 型 指标 
DeepLabv3+ (Xception) DeepLabv3+ (MobileNet) SegNet (ResNet50) U-Net (MobileNet) PSPNet (MobileNet) 
农用 地 precision/% 91.25 89.02 89.46 91.68 89.27 
recall/% 91.16 88.82 88.66 89.31 86.85 
Fl1 值 0.91 0.88 0.89 0.90 0.87 
建筑 用 地 precision/% 81.12 66.73 78.65 80.79 75.83 
recall/% 70.36 47.95 63.46 69.07 71.68 
Fifa 0.74 0.53 0.69 0.74 0.73 
水 体 precision/% 99.24 95.58 97.80 97.76 96.31 
recall/% 90.19 96.37 95.37 95.21 96.77 
FUE 0.94 0.96 0.97 0.96 0.97 
荒漠 precision/% 81.72 80.17 77.16 77.02 79.55 
recall/% 77.64 67.50 74.66 78.30 70.27 
FUE 0.78 0.70 0.75 0.76 0.73 
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#5 非 数 据 集 时 期 影像 的 精度 验证 于 提取 细节 处 的 地 物 信息 :一 是 对 于 小 面积 农用 地 

Tab.S Accuracy verification of images in non-dataset 提取 更 准确 ， 如 景 (BI 21 的 A.B.C 处 ; ; 二 是 对 于 小 
pe 尺寸 细 长 建筑 用 地 的 提取 地 物 网 络 更 加 连贯 ,如 影 
日 期 (年 -月 ) aes 像 块 22 的 A 处 ,影像 块 23 的 A.B 处 ,以 及 影像 块 25 
Kap 的 B 处 ;三 是 对 于 细碎 无 规律 荒漠 地 的 提取 效果 的 
nid 0:35 边缘 细节 更 精细 , 见 影像 块 22 的 B 处 ,影像 块 24、25 
PANOS oie 0.96 的 A 处 ,与 其 余 模 型 的 分 类 结果 相 比 ,相同 位 置 处 的 

_ SPAR 0 0% _ 提取 效果 更 准确 精细 , 误 分 情况 更 少 ,边缘 提取 效 


DeepLabv3+ 
(Xception) 


DeepLabv3+ 


(MobileNet) 


PSPNet 
(MobileNet) 


SegNet 
(ResNet50) 


U-Net 
(MobileNet) 


注 : 各 影像 块 中 A、B、C、D\E\F 为 提取 结果 中 具有 典型 示范 效果 的 局 部 区 域 。 
图 4 各 个 网 络 提取 对 比 结 


Fig.4 Comparison results of each network extracts 
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这 与 表 3 和 表 4 的 地 物 提取 结果 相符 。 

其 余 4 种 网 络 模型 分 别 在 边缘 提取 细小 地 物 、 
地 物 网 络 的 连贯 性 等 方面 表现 不 佳 ,尤其 在 细小 地 
物 提取 方面 ,存在 破坏 地 物 网 络 连贯 性 的 现象 。U- 
Net( MobileNet) 的 编 解 码 架 构 能 够 融合 不 同 层级 特 
征 信 息 ,在 固定 结构 的 医学 影像 上 表现 突出 ,但 是 
在 环境 复杂 背景 的 遥感 影像 上 的 提取 不 精细 ,缺乏 
多 尺度 的 上 下 文 信息 ;SegNet(ResNet50) 边 缘 细 市 
提取 效果 不 理想 因素 包括 对 特征 图 进行 去 池 化 操 
作 时 会 忽略 邻近 信息 ,导致 边缘 轮 廊 模糊 。PSPNet 
(MobileNet) 模 型 没有 结合 丰富 的 空间 信息 进行 解 
码 , 选 用 的 骨干 网 络 MobileNet 的 特征 提取 能 力 不 
强 ,存在 大 面积 椒盐 噪声 ,提取 效果 不 佳 并 存在 误 
分 现象 。 对 比 基 于 Xception 的 DeepLabv3+ 模 型 和 
基于 MobileNet 的 DeepLabv3+ 模 型 总 体 分 类 效果 ， 
验证 了 在 均 采 用 深度 可 分 离 卷 积 结构 情况 下 ,Xcep- 
tion 旨 在 提高 模型 的 性 能 ,而 MobileNet 借 其 减少 参 
数量 提高 计算 速率 ,但 特征 提取 效果 没有 很 大 
提升 。 

3.2 动态 变化 监测 分 析 

3.2.1 总 体 各 地 物 动态 变化 分 析 对 1998 一 2020 
年 22 景 影像 数据 (不 含 2012 年 数据 ) 的 分 类 结果 进 
行 统计 ,总 体 上 ,1998 一 2020 年 农用 地 总 面积 增加 
12.68% ; 荒漠 总 面积 减少 15.00% ;建筑 用 地 面积 总 
体 增加 2.53% ;水 体 总 面积 在 30~40 km’ 范围 上 下 
浮动 。 

FEB PS BK 1998 一 2020 年 地 物 结 构 转 变 明 显 
(图 $) ,尤其 是 农用 地 和 荒漠 2 种 地 物 的 变化 ,建筑 
用 地 和 水 体面 积 均 以 不 同 程度 变化 。 农 用 地 面积 
整体 在 起 伏 间 呈 显 著 上 升 趋势 ,其 中 2013 年 达 最 大 
面积 1059.85 km’; 车 江面 积 呈 明显 下 降 趋势 ,其 中 
2014 4 sie W M AR Ae (ER 308.01 km ,体现 “ 沙 退 人 
进 ” 的 地 物 类 型 变化 规律 ;水 体面 积 结构 没有 明显 
变化 ,结合 分 类 结果 及 实地 考察 ,研究 区 水 体面 积 
直接 影响 因子 为 降水 量 , 水 库 较 浅 的 部 分 在 降水 量 
少 的 时 期 容易 向 荒漠 转移 ,其 中 2017 年 水 体面 积 最 
高 为 31.74 km ,而 2017 年 降水 量 比 历年 低 30%~ 
60%' 引 , 短 时 期 内 自然 环境 因素 的 变化 对 于 其 他 土 
地 覆盖 变化 的 影响 不 显著 ;建筑 用 地 在 2005 年 前 增 
长 趋势 颇 为 缓慢 ,2005 一 2008 年 出 现 较 大 浮动 ， 
2008 年 建筑 用 地 面积 达到 最 小 值 27.95 kim’, Z Ja E 
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图 5 1998 一 2020 年 各 地 物 面积 变化 趋势 
Fig. 5 Change trend of various features area 
from 1998 to 2020 


筑 用 地 面积 以 稳步 趋势 上 升 。 整 体 变化 情况 与 人 
口 增长 .经济 发 展 的 社会 背景 相符 。 
3.2.2 不 同时 段 地 物 结构 变化 分 析 遥感 影像 分 类 
后 ,为 呈现 原本 分 类 效果 ,不 进行 任何 后 处 理 , 将 影 
像 块 反 癌 拼接 制作 地 物 要 素 提 取 专 题 图 。 通 过 
1998 ,2008 .2020 年 专题 图 可 视 化 3 个 研究 时 期 的 各 
类 地 物 面积 变化 情况 (图 6)。4 种 地 物 总 体 分 布 情 
况 : 农 用 地 面积 占 比 最 大 ,荒漠 面积 占 比 次 之 ,主要 
分 布 在 中 北部 区 域 , 水 体 主 要 分 布 在 研究 区 南部 ， 
建筑 用 地 依托 于 农用 地 和 水 域 附近 ,散布 于 整个 研 
究 区 。23 a 以 来 ,研究 区 以 北 荒漠 开明 成 为 农用 地 
等 “ 沙 退 人 进 ” 演 变 过 程 较为 明显 ,由 稀 玲 分 布 的 农 
用 地 逐渐 问 外 扩展 , 变 成 连贯 农用 地 结构 。 水 体面 
积 有 所 浮动 ,但 整体 水 体 区 域 固定 。 建 筑 用 地 呈 扩 
张 趋势 ,以 1998 年 的 建筑 用 地 为 中 心 向 外 扩建 。 
3.2.3 转移 纸 阵 和 动态 度 分 析 借助 1998 一 2020 
年 时 间 序 列 的 土地 利用 转移 矩阵 ( 表 6) ,反映 了 长 
时 间 序 列 内 不 同 地 物 类 型 定量 转化 关系 。 

建筑 用 地 主要 是 向 农用 地 转移 ,转移 比例 达 
97% ,向 水 体 和 荒漠 转移 量 较 小 ;水 体 结构 变化 不 明 
显 , 水 体 主要 是 向 荒漠 转移 ;荒漠 向 农用 地 转移 , 转 
移 比 例 为 95.12% ,其 次 是 向 建筑 用 地 转移 ,向 水 体 
转移 量 较 少 , 近 23 a 英 索 湾 星 区 土地 结构 转变 主要 
方向 是 荒漠 转变 为 绿洲 ,表明 荒漠 化 治理 的 有 效 
性 ;农用 地 主要 是 向 建筑 用 地 转移 ,转移 面积 占 总 
体 转移 面积 的 65.32% ,其 次 是 向 荒漠 转移 ,转移 比 
例 为 25.58% ,农用 地 向 水 体 的 转移 量 较 小 , 且 农 用 
地 疝 建设 用 地 转 入 面积 大 于 转 出 面积 。 
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不 盼 丽 等 ,基于 深度 学 习 的 寒 旱 区 多 时 序 影像 土地 利用 及 变化 监测 一 一 以 新 疆 莫 索 湾 晨 区 为 例 


(a) 1998 年 


地 物 类 型 


(b) 2008 年 
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(c) 2020 年 
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Fig. 6 Comparison of thematic maps of feature extraction in Mosuowan reclamation area in 1998, 2008 and 2020 
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Tab. 6 Transfer matrix of land types in Mosuowan 


reclamation area from 1998 to 2020 /km 
2020 年 1998 年 地 物 类 型 
地 物 类 型 农用 地 ”建筑 用 地 ”水体 i 
农用 地 778.40 11.32 1.02 222.39 
建筑 用 地 37.05 11.61 0.19 11.05 
水 体 0.86 0.01 34.00 0.40 
荒漠 12.81 0.34 3.09 326.52 


漠 ,动态 度 绝对 值 为 1.76% ,以 9.46 km-a EK F 
降 ;总 体 水 体面 积 动态 变化 幅度 不 大 。 在 1998 一 
2008 年 与 2008 一 2020 年 这 2 个 阶段 中 ,农用 地 变化 
程度 在 前 一 个 阶段 增长 比 后 一 个 阶段 剧烈 ,建筑 用 
地 与 之 相反 ,荒漠 面积 的 降低 速率 逐渐 下 降 。 莫 索 
湾 星 区 的 土地 覆盖 及 其 变化 现状 与 国家 宏观 发 展 
导向 息息相关 ,历经 1998 一 2000 年 的 农业 资源 发 展 
阶段 ,2000 一 2010 年 的 西部 大 开发 阶段 ,2010 至 今 


根据 土地 利用 转移 和 矩阵 及 公式 (9) 得 到 莫 索 湾 
垦区 1998—2020 年 地 物 类 型 变化 信息 ( 表 7)。 近 
23 a, 建 筑 用 地 的 地 物 动态 度 最 高 (7.17%) ,以 1.60 
km-a 的 平均 速度 增长 ;农用 地 面积 增长 量 最 大 ， 
面积 平均 增长 8.05 km .a ;面积 下 降 率 最 高 的 是 荒 


为 对 口 援 疆 重 要 阶段 “。 据 兵团 统计 年 鉴 ,1998 一 
2018 年 城镇 化 率 从 56.2% 增 长 为 83.18%, 人 均 CDP 
增长 约 11 倍 强 。2010 年 以 后 ,在 区 域 统筹 实现 兵团 
城乡 一 体 化 发 展 目标 下 , 莫 索 湾 垦 区 加 快 城镇 化 发 
展 步伐 ,新 疆 天 业 .十 户 滩 新 材料 工业 园区 等 工业 
也 迅速 崛起 ,反映 了 经 济 发 展 驱 动 城镇 化 建设 的 现 


表 7 1998 一 2020 年 莫 索 湾 垦区 地 物 类 型 变化 信息 


Tab. 7 Change information of surface feature types in Mosuowan reclamation area from 1998 to 2020 


面积 变化 /km? 单一 土地 利用 动态 度 /% 
地 物 类 型 一 - 
1998—2008 年 2008—2020 4 1998—2020 4E 1998—2008 4E 2008—2020 4E 1998—2020 4 
农用 地 +162.43 +22.77 +185.20 +1.96 +0.23 +1.01 
建筑 用 地 +4.65 +32.12 +36.77 +2.00 +11.49 +7.17 
水 体 -5.49 +2.49 -3.00 -1.43 +0.76 -0.36 
荒漠 -160.86 -56.83 -217.69 -2.87 -1.42 -1.76 


TE: + 表示 增加 ; -表示 减少 。 
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状 , 这 种 地 物 类 型 转变 情况 符合 城镇 扩张 和 经 济 发 
展 的 时 代 背 景 。 


4 结论 


本 文 将 中 分 辨 率 卫星 遥感 技术 与 深度 学 习 算 
法 结合 ,实现 典型 地 物 分 类 ,在 地 物 分 类 结果 的 基 
础 上 借助 土地 利用 转移 矩阵 和 土地 利用 动态 度 对 
土地 覆盖 变化 进行 定量 分 析 ,探索 近 23 a ERS E 
区 的 土地 结构 变化 规律 ,主要 结论 如 下 : 

(1) 相 较 于 传统 机 器 学 习 方法 ,深度 学 习 技 术 
对 地 物 信息 提取 的 地 物 分 割 吻合 度 较 高 并 且 自 动 
化 程度 较 高 ,其 中 DeepLabv3+(Xception) 模 型 提取 
效果 最 优 ,本 文 所 使 用 的 MIoU 较 其 他 模型 提升 了 
0.03~0.39 ; 该 模型 在 多 尺度 地 物 分 制 上 具有 独特 的 
优势 ,在 大 尺度 地 物 (农用 地 Fl 值 为 0.91) .小 尺度 
地 物 (建筑 用 地 Fl 值 为 0.74) 均 能 获得 较 优 的 结 
为 深度 学 习 技术 应 用 于 地 物 分 类 提供 参考 与 借鉴 。 

(2) 总 体 上 1998—2020 EER YS BS X eh 
盖 及 变化 程度 在 不 同 地 物 上 存在 明显 差异 ,具体 呈 
现 如 下 特征 : 近 23 a 农用 地 面积 增加 12.68% ,主要 
增加 在 研究 区 东北 部 绿洲 外 围 ; 荒 漠 面 积 减 少 
15% ,荒漠 向 农用 地 转化 比率 高 达 95.10% ;建筑 用 
地 面积 增加 61.21% ,建筑 用 地 的 动态 度 最 大 
(7.17% ) ,2008 年 后 增 速 明显 ,主要 在 城镇 外 围 的 农 
用 地 (转移 面积 占 比 65.32% ) 和 若 漠 (转移 面积 占 比 
25.58% ) 上 逐步 扩张 ;水 体面 积 在 32.12~39.80 km 
之 间 小 幅 波动 ,主要 影响 因子 是 降水 量 。 安 观 政 
策 .社会 经 济 发 展 和 自然 条 件 因 素 等 共同 决定 莫 索 
湾 明 区 土地 利用 及 覆盖 变化 状况 。 

本 文 挖 气 了 深度 学 习 算 法 在 地 物 提取 分 类 及 
动态 变化 监测 方面 的 潜力 和 优势 ,今后 的 研究 中 ， 
重点 研究 如 何 利用 在 网 络 中 使 用 遥感 影像 的 时 间 
分 状 率 信息 ,并 尝试 在 保证 模型 准确 率 的 同时 提升 
算法 效率 ,进而 实现 对 寒 旱 区 空间 资源 分 布 的 高 巡 
准确 监测 ,以 期 为 国土 管理 部 门 等 提供 可 靠 的 宏观 
决策 数据 。 
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Abstract: In this study, we conducted researchon the extraction of ground elements and the dynamic monitoring 
of land cover change in cold and arid areas with fragile ecological environments. The Mosuowan reclamation 
area in Xinjiang, China was selected as the target area, and the Landsat series satellite images from 2015 to 2019 
were cut into subimages with sizes of 416x416 px. Specifically, the training set, validation set and testset are 306 
subimages, 204 subimages and 120 subimages respectively. Three traditional methods were then evaluated: MLC, 
SVM, and RF. Five semantic segmentation models were also evaluated: DeepLabv3 + (Xception), DeepLabv3 + 
(MobileNet), SegNet (ResNet50), U-Net (MobileNet), and PSPNet (MobileNet). In the evaluation experiments, 
DeepLabv3+ (Xception) was found to achieve the optimal segmentation effect and multiscale feature fusion using 
fewer parameters. The overall accuracy, Kappa coefficient, precision, recall, Fl-score, and MIoU were 96.06%, 
0.96, 87.69% , 83.78% , 0.86, and 0.77, respectively. The MIoU of the DeepLabv3 + (Xception) model was 
significantly better than those of the other four models, improving from 0.03 to 0.39. On the basis of the land use 
classification results of the long-term timeseries remote sensing data from 1998 to 2020, we analyze spatial 
structure change in land use and the associated driving factors. Over the past 23 years, the total areas of desert, 
agricultural land, and construction land have been reduced by 15.00%, 12.68%, 2.53%, respectively. At the same 
time, the amount of water area has remained relatively stable. The overall transformation direction of land use is 
from desert to agricultural land and then from agricultural land to construction land. It can be seen from the 
results that the desertification control was effective, and urbanization rapidly developed. Consequently, this study 
can provide a reference for the application of deep learning in the field of medium-resolution remote sensing 
images, which can be used to realize the dynamic monitoring of land use and change. 

Key words: remote sensing image; long time-series; deep learning; DeepLabv3+; feature classification; dy- 
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